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SUMMARY 


Local heat-transfer coefficients, temperature recovery factors, and 
pressure distributions were measured on a circular cylinder at a nominal 
Mach number of 3*9 over a range of free-stream Reynolds numbers from 
2.1X10 3 to 6.7X10 3 and yaw angles from 0° to 44°. 

It was found that yawing the cylinder reduced the local heat-transfer 
coefficients, the average heat-transfer coefficients, and the pressure 
drag coefficients over the front side of the cylinder. For example, at 
kk° of yaw the average Nusselt number is reduced by 3^ percent and the 
pressure drag by 60 percent. The amount of reduction may be predicted by 
a theory presented herein. Local temperature recovery factors were also 
reduced by yaw, but the amount of reduction is small compared to the 
reduction in heat-transfer coefficients. 

A comparison of these data with other data obtained under widely 
different conditions of body and stream temperature, Mach number, and 
Reynolds number indicates that these factors have little effect upon the 
dropoff of heat transfer due to yaw. 


INTRODUCTION 


Current interest in the flight of aircraft and missiles at high 
supersonic speeds has brought with it the problem of aerodynamic heating 
of the aircraft skin and structure. One of the parts of the aircraft 
where heating is most severe is the leading edge of wings. If these 
leading edges are sharp and thin, there is little material available to 
absorb or dissipate the heat. Also, uneven heating of sharp leading edges 
may result in high thermal stresses. 
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A method of alleviating this problem is to blunt the leading edges 
of wings, which reduces the local rate of heat input compared to a sharp 
leading edge, and provides additional material at the leading edge which 
gives additional strength and increased thermal capacity. 

Blunting the leading edge of a wing normally imposes a drag penalty; 
however, if the leading edge is swept back, the drag due to the blunted 
leading edge can be materially reduced. This fact is demonstrated in 
reference 1 in which drag of a yawed circular cylinder is measured at a 
Mach number of nearly 7* Another advantage to be gained by sweeping 
the leading edge is that the heat transfer rate to the leading edge is 
reduced below that occurring if the leading edge is normal to the direc- 
tion of flight. This benefit is brought about by a reduction in both 
the heat-transfer coefficient and the temperature recovery factor. The 
fact that yawing a circular cylinder reduces the average heat-transfer 
coefficient has been recognized for years by workers in the field of hot- 
wire anemometry. King, in 1914, measured this effect (see ref. 2). Ref- 
erences 3 and 4 summarize later work in this field. Recently, average 
heat-transfer rates to yawed and unyawed wires have been measured (ref. 5) 
and it was found that the reduction of heat transfer by yawing discovered 
by the workers in the field of hot-wire anemometry persisted at Mach num- 
bers of the order of 10. 

Previous experimental work in this field, for the case of supersonic 
flow over the cylinder, has been limited to measurements of average heat- 
transfer coefficients or average heat-transfer rates over either the front 
half of the cylinder or over the entire cylinder. The general purpose of 
the research described in this paper is to study the effect of yaw upon 
both the local heat transfer and the pressure drag of a circular cylinder 
immersed in a supersonic air stream. 

The experimental portion of the investigation consisted of measuring 
local heat-transfer coefficients, local temperature recovery factors, and 
pressure distribution on a 1-inch-diameter circular cylinder at angles of 
yaw from 0° to 44°. The tests were conducted in the Ames 8-inch low- 
density wind tunnel at a nominal Mach number of 3.9 and over a free-stream 
Reynolds number range of from 2.1x10 s to 6.7X10 3 . 

In addition to the experimental portion of the investigation, a 
theory is derived from which local heat-transfer coefficients and pres- 
sure drag coefficients over the front half of a yawed circular cylinder 
may be predicted to an accuracy sufficient for most engineering purposes. 


NOTATION 


a 


speed of sound, ft /sec 
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"D 


E 


T 

constant in relation = C — 

^ a 

pressure drag coefficient, based on projected area 
specific heat of air at constant pressure, ft-lb/slug, °F 
cylinder diameter, ft 

k T 

constant in relation - — = E — 

T t 

new variable in momentum equation ( 5 ) 

Px=o a x=o 


F(A,M) function of yaw angle defined by 


Pt : 




new variable in energy equation defined by - 


i - i. 


g(m) 


k 

M 

Wu 

P 

Pr 

Pt 

<1 

Rm 


P t 2 &00 1 

function of Mach number defined by — — — 77 

Poo H Fi 

local heat-transfer coefficient, Btu/ft 2 , hr, °F 
total enthalpy, ft-lb/slug 

thermal conductivity of air, Btu/ft 2 , hr, °F/ft 

Mach number, dimensionless 

a oo 

Nusselt number, — , dimensionless 
k t 

pressure, lb/ft 2 

c nM’ / \ 

Prandtl number, (consistent units), dimensionless 

k 

wind-tunnel reservoir pressure, microns of mercury absolute 
heat-transfer rate, Btu/hr 


pJJD 

Reynolds number, 


, dimensionless 
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s 

T 

T r 

U 

u,v,w 

*,y, z 

P 

7 

M- 

P 

v 

T) 


nr 

A 

<P 

o(cp) 


p UD 

Reynolds number, , dimensionless 

^t 2 

surface area, ft 2 

t emp erature , °R 

local recovery temperature, °R 

free-stream velocity ahead of normal shock wave, ft/sec 

velocity components in x, y, and z directions, respectively, 
ft/sec 

coordinates on cylinder, ft 

constant of proportionality between velocity Ux and surface 
coordinate x defined by relation u x = (3x 

ratio of specific heats, dimensionless 

viscosity of air, lb-sec/ft 2 

density of air , slugs/ cu ft 

kinematic viscosity, ft 2 /sec 

/py/ 2 

new variable in momentum equation defined by ( "v J z 

T - T 

temperature recovery factor, =^, dimensionless 

i t " 00 

angle of yaw, deg 

azimuth angle measured from forward stagnation point, degrees or 
radians as noted 

Pi 

function of azimuth angle defined by 

■ p x=o 


Subscripts 


surface of body 

conditions at outer edge of boundary layer 
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2 conditions just downstream of normal shock wave 

t total conditions (i.e., conditions that would exist if the gas 

were brought to rest isentropically) 

av average quantity over front half of cylinder 

oo free-stream conditions ahead of shock wave from cylinder 


ANALYSIS 


Before proceeding with the details of the analysis, the main purposes 
will be outlined. Briefly, it was hoped that the theory would yield, as 
a minimum result, correlation parameters or dimensionless groupings which 
could be used to correlate the experimental data, and, secondly, that the 
functional relationships between the local Nusselt number and these parame- 
ters could be deduced. As an additional objective, it was hoped that the 
theory would provide a means by which Nusselt numbers could be predicted 
at flow conditions different from those at which tests have been conducted. 
The degree to which these objectives have been realized is discussed in 
later sections . 



By means of order of magnitude arguments similar to those used by 
Sears (ref. 6) the momentum , continuity, and energy equations for laminar 
flow over an infinite yawed cylinder can be developed. 1 In a coordinate 
system (see sketch (a)) where x is measured along the surface of the 

•’•Crabtree (ref. 7 ) obtains the same set of equations in a recently 
published work. 
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cylinder from the stagnation point in a direction perpendicular to the 
axis, z measured normal to surface and y measured spanwise, the con- 
servation of momentum in the x direction is given by 



Conservation of momentum in y direction reduces to 


pu 


Sv , T T dv 
^ + pw ^ 



(1) 


(2) 


The equation of continuity of mass is given by 



+ ^- (pw) = 0 

oz 


(3) 


The conservation of total energy for a Prandtl number of 1 is given by 


Si 

pu — + pw 
ox 


Si 

Sz 


s / Si\ 

Sl^s -J 


(b) 


where 


1 = 


u 2 + v c 


+ CpT 


If the flow is assumed incompressible with constant properties and 
the following change of variables is made 

U = Pxf^T}) 
w = - (vp-pj ) 1/2 f (ti) 


u x = 0x 
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the momentum equation in the x direction transforms to (see, e.g., 
ref. 8) 

f^ 2 - ff-qq = 1 + f^qq (5) 

and the energy equation transforms to 

®t|T| + f Sq = 0 < 6) 

where the subscript T| denotes differentiation with respect to the new 
variable q . These transformations are given in detail in reference 8. 

The enthalpy variable g(q) is taken to be a function of q only 
and has the limits g( tj) = 1 at q = oo and g( iq ) = 0 at q = 0. It can 
be shown that 8g(q)/8x is a term of small order in comparison to the 
term 8g(q)/8q if the approximations made to derive equations (l) , (2), 
and (4) from the more general Navier-Stokes and energy equation hold. 

The heat-transfer rate per unit area is given by the solution of 
equations ( 5 ) and (6) and is 

f = k(T t - T s ) (£) (gq) s • (7) 


where 


(g„> 3 - 0.570 

The solution to equations ( 5 ) and (6) were obtained by Pohlhausen in 1921 
(see ref. 8) . 

If the heat-transfer coefficient is defined by 


h = 



( 8 ) 


then 


h = O. 57 O k 



(9) 


The above equation ( 9 ) gives no hint as to where in the boundary 
layer the kinematic viscosity and thermal conductivity should be evaluated; 
however, there are various pieces of evidence to guide the choice. One 
method, widely used (see, e.g., ref. 9 ), is to assume a linear relation- 
ship between viscosity and temperature and between thermal conductivity 
and temperature , and to evaluate the pressure at the outer edge of the 
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boundary layer (since 8 p/ 8 z has been assumed equal to 0 in the boundary 
layer) . The constant 3 was evaluated from pressure distribution data 
and will be discussed in a later section. Cohen and Reshotko (ref. 10) 
discuss the effect of Prandtl number on the heat-transfer coefficient. 

They found that the factor Pr 0 * 4 multiplied times the heat-transfer 
coefficient obtained from the analysis where the Prandtl number is assumed 
equal to 1 accounts for this effect. This factor is included in the fol- 

£L 

lowing equation (10a). The above assumptions yield, if 3 = 2.13 ■ X ~° 

D 


h = 0 . 588 ■= Pr° ' 4 k-t 
•J C 



/ y /a /p x = 0 y /2 /a x=0 \^2 

v Pt J v h ) 


(lOa) 


where 


Pi 

Px=o 


1 



(10b) 


Px=Q 

Pt 2 


cos 2 A + 



7 - 1 
7 + 1 / 


1 

7-1 


1 


M 2 


/ 

|7-i 


sin 2 A 


(10c) 


a x=o 

at 



7-1 \ 1/2 

— 2 M 2 sin 2 A\ 



M 2 


(lOd) 


From equation 
coefficient at the 


(10), the effect of yaw angle A upon the heat-transfer 
forward stagnation point is given by 


h A = fPx = oV /2 /a x=0 \ l/2 

^A=o \Pt 2 J V a t / 


( 11 ) 


Also from equation (10) the effect of azimuth angle cp upon the ratio 
of the local heat-transfer coefficients is given by 
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DESCRIPTION OF EQUIPMENT AND EXPERIMENTAL METHOD 

Wind Tunnel 


The tests were conducted in the Ames 8-inch low-density wind tunnel. 
This wind tunnel is an open- jet nonreturn type tunnel. Air was used as 
the test gas. The 8-inch tunnel is a scaled up version of the low- 
density w in d tunnel described in reference 11. A five- stage set of steam 
ejectors is used to produce the main flow. The axisymmetric nozzle was 
designed by the method described in reference 12 , with the additional 
feature of boundary— layer removal. The nozzle was constructed of shim 
stock of varying thickness and alternate shims were removed to permit 
boundary- layer removal as described in reference 13. The design Mach 
number was 4 through the stream-static pressure range of 100 to 400 
microns of mercury. The boundary layer is removed by a set of steam 
ejectors operating in parallel with the main drive ejector set. The 
physical arrangement of the nozzle and test section is shown in figure 1. 

Preliminary surveys of the nozzle indicated that no strong shock 
waves were present in the nozzle when the expansion ratio across it was 
properly set and controlled. The air stream was surveyed with an open- 
end .impact pressure probe. Surveys were made in a plane normal to the 
stream direction l-l/4 inches downstream of the nozzle exit. Surveys 
were also made in the nozzle along the stream center line. The static 
pressure of the stream was obtained by measuring the nozzle wall pressure 
at a point 2 inches upstream of the exit plane of the nozzle. This method 
of obtaining stream static pressure has been described in reference 12. 

A typical Mach number distribution obtained from these measured quan- 
tities is shown in figure 2. The Mach number was calculated in two ways, 
(l) from measured impact pressure and static (wall tap) pressure together 
with the assumption of a normal shock wave standing ahead of the impact 
pressure probe and isentropic deceleration of the flow behind the shock 
wave (circular points) and (2) from measured impact pressure and upstream 
reservoir pressure (total head) using the assumption that the flow through 
the nozzle was isentropic (shown by the square symbols). Good agreement 
was obtained between the two methods of obtaining Mach number over the 
range of pressure levels used in the investigation. Therefore the assump- 
tion that the flow through the nozzle was isentropic (in the test region) 
appears to be reasonable. 

Table I presents the actual usable stream diameter and Mach number 
obtained for various test section static pressures. 
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Model 


The heat-transfer model was a 6-inch-long cylindrical copper shell 
of 1-inch outside diameter and l/^- inch- thick wall (see fig. 3). A copper 
plug, l/8-inch diameter and l/8-inch long, was inserted into a hole in 
this shell, with a l/64-inch air gap between the plug and the shell. The 
surface of the plug was machined to the contour of the cylinder. The body 
(or shell) was instrumented with an electrical heater at each end spaced 
2-l/2 inches from the plug, a thermocouple embedded in the shell under 
each end heater, and a thermocouple in the shell near the plug. An elec- 
trical heating coil was wound on the plug, and a differential thermocouple 
mounted between the plug and the shell. This differential thermocouple 
was used to indicate the temperature difference between the plug and the 
shell. Mechanical means located outside of the stream were provided in 
the mounting to permit rotation and yaw of the cylinder which completely 
spanned the stream. 

A plastic film, 0.00025-inch thick, was wrapped around the cylinder 
to seal the air gap between the plug and the body shell. The air gap was 
then vented to the hollow portion of the cylinder and thence to the con- 
stant static pressure of the test section. Thus the heat loss due to 
conduction through the air gap is reduced as much as possible because of 
the presence of a quiescent layer of low pressure air around the test plug. 

The pressure model was constructed from a 1-inch diameter cylindrical 
shell. A pressure tap of 0.035-inch diameter was located in the center 
of the cylinder, the pressure at the tap was measured by an oil-filled 
U-tube manometer. Pressures were measured for various azimuth positions 
around the cylinder at 0°, 30° , ^5°, an( 3- 60° of yaw. 


Test Method 


The heat transferred from the surface of the plug to the air stream 
was determined as a function of the difference between the plug tempera- 
ture and the stream stagnation temperature. A test point was obtained by 
heating the cylinder and the test plug to the same constant temperature, 
and measuring the plug heater current for this steady-state condition. 

A series of tests were made with no air flow through the tunnel at approxi- 
mately 0.1 micron and at pressures of 100 to 300 microns, to obtain the 
radiation and conduction loss. The variation of these losses with pres- 
sure was within the experimental scatter of the tare data. These heat 
losses were then treated as a tare loss to be subtracted from the gross 
heat input to the plug obtained in the tests. The magnitude of the tares 
was found to be approximately 10 percent of total heat input at highest 
rates, and approximately 60 percent of total heat input at the lowest 
rates present on the back side of the cylinder. At a given orientation 
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of the local test spot, net heat input to the plug was obtained at a 
series of plug temperature levels, ranging from 20° to 50° Fahrenheit 
above stagnation temperature. This net heat into the plug, which is 
then the heat transported from the plug to the stream, was plotted as a 
function of the difference between plug temperature and stagnation tem- 
perature. The slope of this curve is proportional to the product of the 
heat-transfer coefficient and the test area which was taken slightly 
larger than the plug area as explained in Appendix A. Extrapolation of 
the curve to zero heat transported gives an intercept which is the dif- 
ference between recovery temperature and stagnation temperature. A 
typical test curve is shown in figure 4 . Similar experimental curves 
were obtained at azimuth angles of 0° to 90°, at yaw angles of 0°, 30°, 
and 44 °, and for stagnation temperature of 5^0° R. At zero angle of yaw, 
these curves were obtained up to azimuth angles of l80° at one test 
condition. 

Tests were performed on a different body to determine the effect of 
the thickness of the Mylar film covering the plug. Both tare tests and 
heat-transfer tests were made using two different thicknesses of Mylar 
film. The effect of the additional layer of film on the results was 
within the scatter of the data. 


RESULTS AND DISCUSSION 
Experimental Results 


The experimentally determined local heat-transfer coefficients plot- 
ted versus azimuth angle for the case of the cylinder normal to the stream 
are shown in figure 5. The solid lines are faired through the experimental 
points obtained over a range of free- stream Reynolds number. It can be 
seen that heat-transfer coefficients decrease monotonically back to an 
azimuth angle of 90° • Over the rear portion of the cylinder, the heat- 
transfer coefficients are very low compared to the value at the stagnation 
point; the average value being only about 11 percent of its value at the 
stagnation point. 2 It can be seen that lowering the Reynolds number of 
the flow decreases the local heat-transfer coefficients over the front 
half of the cylinder as is also the case in subsonic flow. 

In order to calculate the actual local heat-transfer rates from the 
cylinder, the local recovery temperature must be known. Local free- stream 
temperature recovery factors are shown in figure 6 for the same conditions 

2 The accuracy of the measurement on the rear portion is reduced due tc 
low heat-transfer coefficients and relatively high tares (approximately 
60 percent of total heat input). It may be of interest to point out here 
that the pressures measured on the back side of the cylinder were very 
low, as may be seen from the data tabulated in table II. 
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of flow given in figure 5 . It can be seen that the temperature recovery- 
factor decreases from a value of unity at the forward stagnation point to 
a value of O .67 at an azimuth angle of 120° and then increases toward 
unity as the rearward stagnation point is approached. 3 Also changing the 
Reynolds number of the flow did not appear to materially alter the varia- 
tion of local temperature recovery factor over the front half of the 
cylinder . 

Other investigators have measured local temperature recovery factors 
on cylinders normal to the stream and the results of these tests are sum- 
marized in reference l4. Up to an azimuth angle of 60 the data of this 
reference agreed well with the results of the present tests, as can be 
seen in figure 6. At azimuth angles between 60° and I 5 O 0 the present 
tests give results which are much lower than those of reference 14. For 
azimuth angles between I 5 O 0 and 180°, the results of the present tests 
are higher than those of reference 14. In the experiments described in 
this reference, sharp changes in recovery temperature would tend to average 
out due to heat conduction in the models . 

This dropping-off of recovery factor with azimuth angle tends to 
make the front portion of the cylinder even more controlling of the heat 
rates than would be indicated by the ratios of heat-transfer coefficients 
at the 90° point to those at the forward stagnation point. A statement 
of the amount of heat transferred from the front half compared to that 
transferred by the back half is difficult if not impossible to make unless 
the stagnation temperature of the flow and wall temperatures of the body 
are specified, as the heat-transfer rates depend upon the heat-transfer 
coefficient and the driving temperature potential for all cases where the 
wall temperature is not very small compared to the stream stagnation 
temperature . 

The effect of sweep or yaw angle upon local heat-transfer coefficients 
is shown parametrically in figure 7 wherein the local heat-transfer coef- 
ficient is shown versus azimuth angle for a free- stream Reynolds number of 
6.7X10 3 . Three angles of sweep are shown, 0°, 3°°, and ^°, and it can 
be seen that yawing the cylinder reduces the local heat-transfer coeffi- 
cients at all azimuth angles up to 90°. It is interesting to note . that 
yawing the cylinder reduces the heat-transfer coefficient at any given 
azimuth angle by approximately the same percentage. 


As companion information with figure 7, the local temperature recovery 
factors at the aforementioned angles of yaw are shown in figure 8 as a 
function of the azimuth angle. It can be seen from figure 8 that yawing 
the cylinder reduces the temperature recovery factors. The reduction in 
recovery factor is, however, small compared to the reduction in heat- 
transfer coefficient produced by yawing the cylinder. The heat-transfer 
data are tabulated in table III . 


3 See footnote 2, p . 11. 
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Comparison of Experimental Results and Analysis 


Pressure distribution and drag .- Three major assumptions had to he 
made in the analysis in order to simplify the basic differential equations 
governing the flow sufficiently to allow a solution. These assumptions 
were that the Prandtl number was equal to 1, that the flow was incompressible 
and properties were constant, and that the x component of the external 
velocity over the front half of the cylinder could be expressed as u x = (3x. 

The net result of the assumption of Pr = 1 is that the analysis 
yields a recovery factor of 1 or that the recovery temperature of the 
cylinder is constant and equal to the stream stagnation temperature. That 
this is not the case can be seen from the data in figures 6 and 8. A simi- 
lar difficulty arises when this assumption is made in analysis of flow over 
flat plates. Experience has shown, for the case of flow over flat plates, 
that the assumption of Pr = 1 causes the calculated Nusselt number to be 
higher than the experimental value by a constant factor equal to Pr 1/3 . 
When the theoretical value of the Nusselt number, obtained by assuming 
Pr = 1, is used to calculate the actual heat -transfer rate from a flat 
plate it must be multiplied by Pr 1 ^ 3 and the experimental value of the 
recovery temperature must be used in the temperature potential in order to 
obtain results that are in agreement with theory. Cohen and Reshotko 
(ref. 10) discuss the factor Pr 0,4 used to correct the theoretical Nusselt 
number, obtained by assuming Pr = 1, for unyawed two-dimensional bodies. 
The assumption is made here that this factor applies to the yawed cylinder 
as well. 

The assumption that the flow was incompressible and that properties 
were constant is probably the weakest assumption made in the analysis. 
However, in the application of the analysis the viscosity is allowed to 
vary linearly with temperature, and pressure is evaluated at the point on 
the surface being considered. Comparison of the results of the analysis 
with experiment will be made in a later section to check the validity of 
this assumption. 

The assumption that the x component of the external velocity over 
the front half of the cylinder is a linear function of x was checked in 
the following way. The velocity over the cylinder in the x direction 
was calculated from measured pressure distribution using Bernoulli’s equa- 
tion for a compressible gas and assuming the fluid velocity was zero at 
the stagnation point. When this was done it was found that if the con- 
stant p was set equal to 2.13 a x=0 /D the velocity over the cylinder, 
yawed or unyawed, could he calculated with good accuracy . By substituting 
this expression into Bernoulli’s equation, the pressure distribution over 
the front half of the cylinder could be calculated. Figure 9 shows the 
ratio of the pressure at any azimuth angle to that at the forward stagna- 
tion point plotted versus azimuth angle. It can be seen from figure 9 
that, for the case of a cylinder normal to the stream, variation in 
Reynolds number from 6.7X10 3 to 1 . 4x10 5 and variation in Mach number 
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from 2.48 to 6.86 does not appreciably alter the pressure ratio distri- 
bution over the front half of the cylinder. The solid line was obtained 
by substituting u x = px into Bernoulli's equation and is given by 


P 


Pi 

x=o 


4 



(13) 


Figure 10 shows the same parameters as figure 9; however, the test points 
were obtained at yaw angles from 0° to 60°. The solid line is again the 
curve calculated from equation (13). Although yawing the cylinder does 
change the pressure over it, from figure 10 it can he seen that the pres- 
sure ratio variation is not changed for yaw angle of 0° to 30° tn the 
present tesbs and 0° to 60° in the tests of reference 1. The pressure 
distribution over the cylinder measured at 45° and 60° of yaw in the pres- 
ent tests departed from that reported in reference 1 for azimuth angles 
greater than about 45°. It is suspected that this departure is brought 
about by the fact that the flow over the cylinder was becoming three- 
dimensional due to the 1-inch-diameter model in the 3-inch-diameter stream. 
The conclusion can then be drawn that over the range of variables investi- 
gated (R = 6.7X10 3 to 1 . 4xl0 5 , M = 2.48 to 6.86) that the pressure ratio 
distribution is a unique function of the azimuth angle for yaw angles of 
0° to 30° for the present tests and 0° to 60° for the tests of reference 1. 
Thus, the assumption of u x = px appears to be a reasonable one. 


One other assumption must be investigated before the results of the 
analysis are compared with the experimentally determined heat-transfer 
results, namely, that the pressure at x = 0 on a yawed cylinder may be 
computed by Rayliegh 1 s equation using the component of the Mach number 
normal to the cylinder. Figure 11 shows the ratio of the pressure at 
x = 0 to the stream impact pressure plotted versus the yaw angle of the 
cylinder. The curves were calculated for three Mach numbers using the 
above-mentioned assumption and the test points are from the present tests 
and from reference 1. The good agreement between the curves calculated 
by equation (10c) and the test points indicate that this assumption is 
also a reasonable one. 


The pressure distribution over the front side of a yawed or unyawed 
cylinder can now be used to compute the effect of yaw upon the pressure 
drag over the front side of the cylinder. The resultant expression is 


^A 

(C D ) A=o 


Px ° cos A 
P t 2 


(14) 


Figure 12 shows this ratio plotted as a function of the angle of yaw of 
the cylinder for two Mach numbers. The curves are obtained from equa- 
tion (l4) and the test points shown are from the present tests and from 
reference 1. The agreement between equation (l4) and the experimental 
points is good. This figure points up the previously mentioned fact 
that rounding the leading edge of a wing may not lead to severe drag 
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penalties if the wing is swept. For example, at 45° of sweep the drag 
coefficient of the swept leading edge is only about 40 percent of its 
unswept value. Also, the drag coefficient ratio of equation (l4) 
approaches cos 3 A for Mach numbers approaching infinity, which corre- 
sponds to Newtonian flow results. 

Local Nusselt number .- The ratio of the local Nusselt number or 
local heat-transfer coefficient to that at the stagnation point is plot- 
ted versus azimuth angle in figure 13* Included in this figure are data 
for three Reynolds numbers and three angles of yaw, 0°, 30°, and 44°. It 
can be seen that in general the heat-transfer coefficient decreases with 
azimuth angle. Except for the data taken at 44° of yaw, all of the points 
tend to form a single curve. This fact tends to substantiate the result 
of the analysis which shows that this ratio is a function of azimuth angle 
only. The variation of heat-transfer coefficient ratio given by the analy- 
sis is shown as the solid curve and was calculated from the expression 


It can be seen that at the lower azimuth angles the analysis fits the 
data reasonably well, but at the azimuth angles of 60° and 75° the theory 
predicts values larger than those observed experimentally. 

At an azimuth angle of 60° the Mach number at the edge of the boundary 
layer has reached a value of about 1.2 and at the 75° point the Mach num- 
ber is I. 65 . 

In order to determine if compressibility was responsible for the 
dropoff of the measured values of heat transfer below those given by the 
analysis for these azimuth angles, the theory of Cohen and Reshotko, 
reference 15, was compared with the data for the case of zero yaw. This 
theory, which accounts for the effects of compressibility but not for the 
effects of yaw, is shown by the dashed curves. It is apparent that the 
theory fits the data better at these higher azimuth angles than the incom- 
pressible one. However, if average values of the Nusselt number are con- 
sidered, the difference between the compressible and the incompressible 
theory is a constant, and because of the uniqueness of the Mach number dis- 
tribution over the front half of the cylinder, the incompressible theory 
may be used to correlate data over a wide range of Mach numbers. 

The result of the present analysis (eq. (10)) may be written in terms 
of the local Nusselt number, Prandtl number, Reynolds number evaluated 
behind the bow shock wave, a function of the free- stream Mach number, a 
function of the azimuth angle, and a function of both the yaw angle and 
the free-stream Mach number. The local Nusselt number is then given by 



( 15 ) 


Wu local = tt = °* 8 3 2 Pr o - 4 N^N/F(A,M)G(M)0(cp) 


( 16 ) 
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where 


■ (fe) (^if) 

G(M) - Gr) («?) (h) 

®(cp) = 

F x=o 

Equation ( 16 ) is compared with the data on a local basis in figure 14 
where the local Nusselt number is shown plotted versus the parameter 
R 2 F(A,M)g(m)$(cp). The result of the analysis is shown as the solid curve 
and the test points shown are for three Reynolds numbers, angles of yaw 
from 0° to 44°, and azimuth angles from 0° to 60°. The data are correlated 
by the analysis to within a mean deviation of 10 percent. 

Average Nusselt number .- The result of the analysis is compared with 
the experimental data in figure 15 whereon the average Nusselt number for 
the front half of the cylinder is plotted versus the dimensionless parame- 
ter R 2 F(A,M)g(m) . The solid line is the result of the analysis and is 
given by the following expression which was obtained by integrating $(cp) 
over the front half of the cylinder 

Nu av = O .5935 Pr °' 4 n/r^ [F(A,M)G(M )] 1/2 ( 17 ) 

where Nu av and R 2 are evaluated using free- stream density and velocity 
but viscosity and conductivity evaluated at stagnation conditions. 

The experimental points in the figure (solid points) were obtained 
during the present investigation at three Reynolds numbers and at angles 
of yaw of 0 °, 30 °, and 44°. 

The agreement between the calculated and experimental values of aver- 
age Nusselt number over the front half of the cylinder is within ±10 per- 
cent for all points except the value obtained at zero angle of yaw at the 
lowest Reynolds number (r 2 = 6l0) . This point is 15 percent below the 
predicted curve. 

Also shown in figure 15 are average Nusselt numbers for the front 
half of yawed and unyawed cylinders obtained in the 11 - inch wind tunnel 
at Langley Field, reference 16, at a Mach number of 6.9 and at free-stream 
Reynolds numbers of I. 3 XIO 5 and 1.8xl0 5 . The yaw angle was varied from 
0° to 75° in these tests. It can be seen that the results of the analysis 
correlate the Mach number 6.9 data reasonably well up to angles of yaw of 
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60°. At an angle of yaw of 75° "the data of reference 1 6 do not correlate 
well with the result of the analysis. It is stated in reference 16 that 
the data taken at 75° °f yaw may not he reliable due to model limitations. 
This effect is attributed to lack of two-dimensionality at the high yaw 
angle. The data of reference 16 were obtained with heat flow into the 
model at a stagnation temperature of ll40° R and over a range of model 
temperatures from 570° R to 910° R. The present data were obtained with 
heat flow out of a model at a stagnation temperature of 520° R and model 
temperatures of 540° R to 570° R. No effect on the heat-transfer results 
could be detected under these widely different conditions. 

The effect of yaw upon the average Nusselt number over the front 
half of the cylinder can best be shown in the next figure (fig. 16) where 
the ratio of Nusselt number obtained at yaw to that obtained at zero yaw 
is plotted as a function of the angle of yaw. Also shown in this figure 
are data from reference 16. It can be seen that yawing a cylinder reduces 
the average Nusselt number over that obtained at zero yaw. At 3 O 0 °f y aw 
the reduction shown by the present data is approximately l6 percent and 
at 44°, 33 percent of the zero yaw value. The curves shown in the figure 
are the result of the analysis and were calculated from the following 
expression 


(Nu a v) A 

h A=o (Nu av ) A=Q 


[f(a,m) ] 1/2 


(18) 


for three Mach numbers, 4, 7, and 00 . It can be seen that F(A,M) is a 
weak function of the free-stream Mach number at yaw angles less than about 
45 ° but for large yaw angles the theory predicts a sizeable effect of 


Mach number on 


n A 


tyv=o 


The data and the predicted result are in good agree- 


ment up to angles of yaw of 44°. At an angle of yaw of 60° and 75 , the 
data of reference 16 lie above the predicted curve. 


An effort was made in the present investigation to extend the range 
of the tests to an angle of yaw of 60°; however, an examination of the 
pressure distribution over the model at this yaw angle (see fig. 10) dis- 
closed a departure from that obtained at the lower yaw angles. At the 
lower yaw angles, namely 30 ° > an examination of figure 10 reveals that 
the pressure distribution over the cylinder agreed very well with that 
reported in reference 1 where the flow was shown to be two-dimensional. 

It is suspected that the deviation at 44° and 60° yaw angles was due to 
the flow over the cylinder becoming three-dimensional because of the 
relatively large model ( 1-inch diameter) in the 3 -inch- diameter stream. 

Heat-transfer results obtained at 60° of yaw also exhibited large 
scatter (about 37-percent maximum spread) and an examination of the model 
revealed that small air bubbles were present between the measuring plug 
and the plastic film. Also electrical shorts between the cylinder test 
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body and the plug heater wires developed about this time. For these 
reasons the data obtained at 60° of yaw were considered unreliable and 
are not included. 


Application to Flight Conditions 


The results, obtained during the present investigation, were for the 
case where the body temperature was very nearly the stream stagnation 
temperature. At high Mach number (5 or above) the stagnation temperature 
obtained during flight may be so high that the aircraft or missile must 
be cooled to a temperature much below the stream stagnation temperature. 

It is of interest, then, to compare the results of these tests and this 
analysis with any data which are available that approximate (as far as 
temperature is concerned) flight conditions. The results reported in 
reference 5 were obtained at a Mach number of 9.8, a stream stagnation 
temperature of 2200° R and a test body temperature of 520° R. The tests 
were conducted on small wires (0.003- to 0.020- inch diameter) at angles 
of yaw up to 70° • 

Even though the tests of reference 5 were conducted at a relatively 
low Reynolds number (315 for the 0.003-inch wire), the data when compared 
with the present tests should indicate in a limited way whether the results 
of the present tests may be applied to the case of a cool body in a hot 
hypersonic air stream. In reference 5 the recovery temperature could not 
be measured; therefore, comparison will have to be made by applying the 
results of the present tests to the specific conditions -under which the 
experiments reported in reference 5 were made. The results reported in 
reference 5 were for heat transfer from the entire cylinder, whereas the 
present tests are for the front side of the cylinder only. During one 
run in the present tests the local heat transfer on the back side of the 
cylinder for the zero yaw condition was obtained. These results are shown 
in figure 7 and it can be seen that the heat-transfer coefficients on the 
back side of the cylinder were low, being only about 10 percent of the 
value at the stagnation point. In the comparison shown in figure 17 it 
was assumed that the heat transfer from the back side of the cylinder 
reported in reference 5 could be neglected. 

Figure 17 shows such a comparison whereon the ratio of total heat- 
transfer rate at angle of yaw to that obtained at zero angle of yaw is 
plotted versus angle of yaw. The circled symbols are from reference 5 
and the square symbols are the results of the present tests applied to 
the above-mentioned stream conditions. The agreement between the two sets 
of data is good, and within the scatter of the data there again appears to 
be no effect of temperature potential upon the dropoff of heat transfer 
with yaw. The solid line is obtained from equation (l8) with the assump- 
tion of constant temperature recovery factor equal to 1. 


NACA RM A55H31 


19 


CONCLUSIONS 


The following general conclusions are drawn from the results of this 
investigation: 

1. Local heat-transfer coefficients, average heat-transfer coeffi- 
cients, and pressure drag coefficients for the front side of a circular 
cylinder are reduced by yawing the cylinder as found by other investi- 
gators, For example, at 44° of yaw the average Nusselt number is reduced 
"by 34 percent and the pressure drag by 60 percent. The amount of the 
reduction may be predicted with sufficient accuracy for most engineering 
purposes by a theory presented herein. 

2. Local temperature recovery factors on the front side of a cylinder 
are reduced by yaw. But this effect is small compared to the reduction in 
heat-transfer coefficients. 

3 . A comparison of these data obtained with body temperature near 
stream stagnation temperature with other data obtained with a varying body 
temperature in a hot hypersonic air stream indicates that these widely 
different temperature conditions have little effect upon the dropoff of 
heat transfer due to yaw. 

4. The heat-transfer coefficients on the back side of a cylinder 
normal to the stream were insignificant compared to those on the front 
side for M = 3*9 and a free- stream Reynolds number of 6.7X10 3 * 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 31 , 1955 
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APPENDIX A 

DETERMINATION OF THE HEAT -TRANSFER AREA 


The net heat was considered to he transferred from the top surface 
of the cylindrical test plug to the air stream. Thus, the area, S, used 
in the following heat-transfer rate equation, was the area of the top of 
the test plug. 


q = hS(T - T r ) (Al) 

However, the film stretched over the model (see fig. 3 insert) does con- 
duct some heat away from the plug. Also, the film receives energy from 
the cylindrical surface of the plug by free molecular conduction through 
the annular air space around the plug. This film acts much as a circular 
fin in dissipating the heat from the test plug and hence it is necessary 
to increase the area to be used in calculating the heat-transfer coeffi- 
cient. The following analysis of the fin effect leads to a determination 
of the correction to the test area. 

The differential equation governing the temperature distribution in 
the circular-film fin may be found by summing the quantities of heat 
transferred by the various means to and from a circular element of the 
fin. Azimuth variations of these quantities around the test plug will be 
considered negligible. The assumption is made that temperature differences 
are small, so that the radiation exchange terms may be written in linear 
form. The width of the annular air space is of the order of a mean-free- 
path length of the gas; thus, it is assumed that the circular element of 
fin gains heat from the plug by free molecular conduction through the 
annular air space. The elemental fin also transfers heat to the stream 
by convection through the flow boundary layer. A further assumption will 
be made that the variation of heat-transfer coefficient and recovery 
temperature is negligible over the area of film considered. On the basis 
of these assumptions and with normal conduction in the film fin, the dif- 
ferential equation is found to be, 
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where 

T temperature along the fin radius 

r radial distance from center of plug 

and the constants are defined by 


. /3MV m * 

V 2-TH 


UcreTo 3 




3^rVm K 
- 2*/jT 


+ 4creT 0 3 



(A3) 


(a4) 


and 

kj thermal conductivity of film, Btu/hr ft 2 °F/ft 
t thickness of film, ft 

N number of molecules per unit volume , l/ft 3 

V m most probable molecular speed, ft/sec 

k Boltzmann constant, 7*2 7*10 -27 Btu/molecule °F 

0 Stefan-Boltzmann constant, 4.8xl0 13 Btu/ft 2 sec °R 4 

e emissivity, dimensionless 

T p temperature of plug, °R 

To temperature of surrounding surfaces , °R 

This differential equation (A2) is a form of Bessel’s equation. The 
solution may be written in terms of modified Bessel functions of zero 
order, first and second kinds, as 

T - — = A 3 Io(r\/B) + A 4 K G (r n/b) (A5) 

B 

where A 3 and A 4 are constants of integration to be determined by the 
following boundary conditions. 
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r = r p , T = T p ' 

r = r B , T = T P I 


r — r — 0 

r - r m, dr - 0 


where 

rp radius of plug, ft 

rp outer radius of annular space around plug, ft 

r m radius at which minimum temperature occurs on fin, ft 


(A 6) 


Since the value of r m is not known, the three boundary conditions (A6) 
determine r m as well as the two constants of integration. 


In the range of interest, the modified Bessel functions in the solu- 
tion (A5) may be replaced by the asymptotic expressions (ref. 17) for 
large values of the argument r n/"b . These expressions are 


lo (r n/"b ) = EMkM 

J 2jtr 7 b 

t 

KoCrjB - ) = I— ■— exp(-r sTb ) 

■7 2 rVB 


(A7) 


We introduce relations (A7) into equation (A5) , apply the boundary con- 
ditions, and evaluate the constants of integration. The radius, r m , at 
which the minimum temperature occurs is found to be very nearly the aver- 
age radius given by 


r m = 


r B + r P 


(A8) 


The temperature distribution in the circular fin is then given by the 
resulting form of equation (A5) as 


fl p 

exp(r'TlB) + 1 

V 2 r m 7 B + l/ 

exp(2r m '/i -t-Jb) 

J r 

exp(r-p\lB) + | 

^2r m v/i-l> 

) exp ( - rp •Tb) 

V2r m VB + l/' 


(A9) 
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The heat removed from the plug by the film is transferred to the 
stream according to the following relation 

dq = h(T - T r )dS (AlO) 

This expression is integrated using the temperature distribution found 
in equation (A9) for T and assuming a constant heat-transfer coefficient. 
The resulting expression is 


q = 2jth 



T r )r dr 


(All) 


If we assume that the amount of heat represented in equation (All) were 
to be transferred at plug temperature from an area given by an equivalent 
radius, r, we have 


q = 2jth(r 2 - rp 2 ) (A12) 

Thus by equating the right-hand sides of equations (All) and (A12) we 
can express the heat-transfer radius r as follows: 



In actual computation of the correction, the emissivity of the filmi 
was taken to be 0.1 (i.e., that of the chrome-plated plug surface). The 
emissivity of the plug with film was found to check closely with the value 
normally taken for polished chrome. The conductivity of the film was taken 
as 0.1 Btu/hr ft 2 °F/ft. This value was obtained from the manufacturer’s 
literature, and was not checked experimentally during these tests. The 
first approximation to h, found by using rp, was used to determine r. 

The correct heat-transfer area is then found using the r computed from 
equation (AI3) . The correction to the area of the plug is approximately 
20 percent for the tests reported herein. 
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TABLE I.- STREAM CONDITIONS 


Static pressure, 
microns Hg abs 

Mach number 

Stream diameter, 
in. 

320 

3.90 

3.6 

180 

3.90 

3-0 

110 

3.75 

3.0 


TABLE II.- PRESSURE SURVEY DATA 


A, 

deg 

<P, 

deg 

Pi. 

Hg abs 

Poo, 

micron Hg abs 

M 

r 2 

R 00 

0 

0 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 

180 

6.54 

6.18 

5.17 

3.80 

2.48 

1.43 

• 75 
.42 
.28 
.24 

.25 

.26 

.26 

.319 

3.94 

1890 

6 . 7 X 10 3 

30 

0 

15 

30 

45 

60 

75 

90 

5.05 

4.77 

3.97 

2.80 

1.85 

1.08 

.60 

.318 

3-90 

1890 

6.7X10 3 

45 

0 

15 

30 

45 

60 

75 

90 

1.05 
1.00 
.85 
.66 
.46 
• 29 
.18 

.10 

3-8 

610 

2 . 1 X 10 3 
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TABLE III.- TEST DATA 


A, 

deg 

deg 

h, 

Btu 

S, 

10- 4 ft 2 

Tr, 

°R 

T t> 

°R 

M 

Poo. 

microns 
Hg abs 

r 2 

Roo 

hr ft 2 °F 

0 

0 

9.05 

1.012 

532 

534 

3.94 

306 

1890 

6.7xio 3 


15 

8.78 

1.013 

532 

534 






30 

7-86 

1.018 

526 

532 






45 

6.46 

1.025 

519 

532 






60 

4.6 2 

1.032 

504 

532 






75 

3.21 

i.o4o 

489 

535 






90 

I.83 

1.059 

459 

537 






105 

1.14 

1.066 

441 

537 






120 

.69 

1.070 

402 

534 






135 

.63 

1.070 

438 

533 






150 

.44 

1.075 

498 

534 






165 

.47 

1.072 

534 

534 






180 

.58 

1.070 

534 

5.34 






0 

Ol — 

1.025 

526 

529 

3.80 

105 

610 

2.lxio 3 


15 

4.63 

1.024 

526 

529 






30 

3.86 

1.030 

525 

535 






45 

3.41 

1.033 

520 

535 






60 

2.52 

1.042 

509 

535 






75 

1.68 

1.052 

472 

529 






90 

1.22 

1.040 

467 

529 






— 

7.77 

1.010 

524 

523 

3.90 

180 

1120 

3 .8xio 3 


15 

7.51 

1.012 

525 

525 






30 

6.50 

1.017 

521 

526 






45 

5.24 

1.024 

506 

525 






60 

3.85 

1.035 

500 

527 






75 

2.57 

1.047 

480 

528 






90 

1.56 

1.058 

457 

528 





30 

0 

TM — 

1.019 

519 

532 

3.91 

300 

O 

ON 

OO 

1 — 1 

6.3XIO 3 


15 

7.59 

1.018 

519 

532 






30 

6.82 

1.023 

511 

527 






45 

5.44 

1.031 

503 

527 






60 

3.43 

1.045 

495 

531 



% 



75 

2.39 

1.053 

473 

533 






90 

1.24 

1.064 

434 

533 






0 

— 03 

1.026 

509 

530 

3.91 

320 

1890 

CO 

O 

1 — 1 

X 

t>- 

• 

xO 


15 

6.12 

1.028 

506 

530 






30 

4.75 

1.036 

489 

522 






45 

3*97 

1.042 

474 

522 






60 

2.75 

1.051 

467 

523 






75 

1.96 

1.058 

453 

525 






90 

1.45 

1.063 

460 

530 









Figure I. - General arrangement of wind-tunnel test section, nozzle, and model. 
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Distance from center of stream, inches 

Figure 2.- Variation of Mach number with distance from center of stream 
for Reynolds number per foot of 8.7XIO 4 at axial distance of I.25 
inches from exit plane of nozzle. 
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End heater 


Heat transferred, q, 10 3 Btu/hr 
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-20 0 20 40 60 80 

Temperature difference, T s -T t ,°F 

Figure 4.- Variation of heat transferred from surface test area to the 
stream with difference between test area temperature and stream stag- 
nation temperature for M = 3-94; Roo = 6. 7X10 3 , cylinder normal to 
stream and test area oriented 15° from the stagnation point . 


Heat-transfer coefficient, 


cn 

3 



Figure 5.- Variation of local heat-transfer coefficient with azimuth angle for a cylinder at zero 

angle of yaw and for three free- stream Reynolds numbers. 
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Temperature recovery factor, rj 
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Figure 6.- Variation of temperature recovery factor with azimuth angle for a cylinder at 

angle of yaw and for three free-stream Reynolds numbers. 


zero 


NACA RM A55H31 


Heat- transfer coefficient, 



Figure 7«~ Variation of local heat-transfer coefficient with azimuth angle for a cylinder at 
various angles of yaw, for a free -stream Reynolds number of 6.7X10 3 . 
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Temperature recovery factor, r) r , 


oo 

O' 



I 

w 

I— 



Azimuth angle, <j>, deg 


Figure 8.- Variation of temperature recovery factor with azimuth angle for a cylinder at various 

angles of yaw for a free-stream Reynolds number of 6.7X10 3 . 
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9.- Variation of the ratio of local to stagnation-point pressure with azimuth angle for a 

cylinder at zero angle of yaw. 
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Pressure ratio, P|/p x = 


OO 

Co 

Azimuth angle, 2x/D, deg 

10 20 30 40 50 60 70 80 90 

|.|, , 1 — , 1 , — 1 >-* \ *-1 


Yaw angle Source 

Q 30 Present tests, M = 3.9 
6 45 



0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 


Azimuth angle, 2x/D, radians 

Figure 10.- Variation of the ratio of local pressure to the pressure at x = 0 with azimuth angle 

for a cylinder at various angles of yaw. 
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Pressure ratio 
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Figure 11.- Variation of ratio of pressure at x = 0 to stream impact pressure with angle of yaw 

of a cylinder. 
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Drag coefficient ratio, 


-r 

o 




Yaw angle, deg 


Figure 12.- Variation of ratio of pressure drag coefficient for front side of yawed to unyawed 

cylinders with angle of yaw. 
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Figure 13.- Theoretical and e^erimental variation of the ratio of the local heat-transfer coef- 
ficient to heat-transfer coefficient at cp = 0 with azimuth angle for various angles of yaw 
and Reynolds numbers. 
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Figure l4.- 
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Comparison of experimental and theoretical variation of local Nusselt number with the 
R 2 F(A,M)G(M)o(cp) for various angles of yaw, azimuth angles, and Reynolds numbers. 
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Figure 15 
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Comparison of experimental and predicted variation of Nusselt number, averaged over 
the front half of the cylinder, with the quantity R^F(A,M)G(M) . 
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Figure 1 6 .- Comparison of experimental and predicted variation of ratio of the average Nusselt 
number for yawed cylinders to average Nusselt number of cylinder at zero yaw. 


NACA RM A55H31 


Ratio of heat-transfer rates, q A /q A = 


z 

> 

o 

> 


r 



Figure 17*- Comparison of experimental and predicted variation of the ratio of total heat trans- 
ferred from front side of cylinder with angle of yaw; M = 9 . 8 , = 2200° R. Data from refer- 

ence 5 include heat transferred from back side of cylinder. 
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